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Abstract 

Paneling  methods  are  approximate  techniques  for  solving  flow 
problems  over  wings  and  bodies.  Vortex  panels  are  used  to  model  flow 
over  wings  and  other  lifting  surfaces.  The  author  develops  a  triangular 
vortex  panel  having  a  vortlclty  distribution  that  can  vary  In  magnitude 
and  direction.  This  panel  Is  used  to  predict  the  pressure  distribution 
on  a  rectangular  and  a  sweptback  wing  in  subsonic  flow.  Lift  distribu¬ 
tions  obtained  compare  favorably  to  Anderson's  solution  and  wind  tunnel 
results  except  near  the  wing  tip.  In  this  region,  the  distribution  will 
spike  before  satisfying  the  Kutta  condition  Imposed  at  the  tip. 

Possible  remedies  for  the  tip  problem  are  discussed.  , 
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DEVELOPMENT  AND  APPLICATION 
OF  A  SUBSONIC  TRIANGULAR 
VORTEX  PANEL 


I.  Introduction 


Background 

Paneling  methods  are  approximate  techniques  for  solving  linearized 
subsonic  and  supersonic  potential  flow  problems  over  wings  and  bodies. 
Panels  In  use  today  Incorporate  singularity  distributions  of  the 
source,  vortex,  and  doublet  types.  Source  panels  are  often  used  to 
model  bodies  **.^d  other  non-lifting  surfaces  and  to  model  thickness 
of  lifting  surfaces.  Vortex  panels  are  used  to  model  either  lifting  or 
non-llftlng  surfaces.  In  a  typical  problem,  the  airplane  Is  represented 
by  a  finite  number  of  panels.  Each  panel  is  a  singularity  distribution 
of  unknown  strength  that  models  some  part  of  the  aerodynamic  surface. 
These  unknown  strengths  are  determined  by  applying  the  flow  tangency 
boundary  conditions  at  or  near  the  aerodynamic  surfaces.  Once  the 
strengths  are  known,  the  perturbation  velocities  can  be  computed.  These 
are  substituted  Into  the  Bernoulli  equation  to  obtain  the  pressure 
distribution  and  the  corresponding  aerodynamic  forces  and  moments. 

Problem  Statement 

In  many  current  paneling  methods,  the  orientation  of  the  vortlclty 
vector  Is  fixed  on  the  panel  (Ref  6  and  Ref  8) .  This  leads  to 
unacceptable  errors  In  some  cases.  The  purpose  of  this  study  is  to 
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derive  a  subsonic  triangular  panel  having  a  vorticity  distribution  that 
can  vary  both  In  magnitude  and  direction.  A  panel  system  Is  assembled 
and  used  to  predict  the  pressure  distribution  on  a  planform.  Since  the 
computations  involved  are  too  lengthy  to  be  performed  by  hand,  a 
computer  code  was  developed  to  apply  the  panels  to  planar  wings. 

Approach 

The  subsonic  triangular  panel  is  derived  in  Section  II.  The 
Biot-Savart  Lav  is  used  to  compute  the  induced  velocity  at  a  point  due 
to  an  assumed  bilinear  vorticity  distribution  on  the  panel.  This 
expression  formulates  the  induced  velocity  in  terms  of  panel  geometry 
and  unknown  corner  point  vorticity. 

Section  III  presents  methodology  for  panel  system  assembly. 
Numbering  Bchemes  are  developed  for  panels  and  the  unknown  corner  point 
vortlclties.  Planform  boundary  conditions  are  applied  which  reduce  the 
number  of  unknowns  and  the  remaining  unknowns  are  solved  for  by 
formulating  a  linear  system  of  equations.  This  system  consists  of 
control  point  equations  (one  per  panel)  and  a  number  of  edge  continuity 
conditions.  The  linearized  form  of  the  flow  tangency  boundary  condi¬ 
tion  is  then  used  to  effect  the  solution.  Once  the  corner  vortlclties 
are  known,  the  vorticity  at  any  point  on  the  planform  can  be  obtained. 
Finally,  induced  velocities  and  corresponding  pressures  can  be  calcu¬ 
lated  from  the  known  vorticity. 

Section  IV  presents  a  computer  code  developed  to  apply  the 
methodology  to  planar  wings.  Brief  descriptions  of  each  subroutine 
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and  a  detailed  Input  description  are  provided.  Appendix  B  is  a 
sample  output  and  Appendix  C  is  the  program  listing. 

Section  V  presents  program  predictions  for  a  rectangular  wing 
which  are  compared  against  Anderson's  solution  (Ref  1:9-16).  Pre¬ 
dictions  are  also  presented  for  a  swept  untapered  wing  which  are 
compared  against  wind  tunnel  tests  (Ref  4:92). 

Section  VI  concludes  the  report  and  makes  recommendations  for  the 
Improvement  of  the  aerodynamic  model. 
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II.  Panel  Derivation 


This  section  presents  the  development  of  the  subsonic  triangular 
panel.  The  goal  is  to  derive  an  expression  for  the  Induced  velocity  at 
an  arbitrarily  chosen  point  In  the  xy  plane  due  to  an  assumed  bilinear 
vorticity  distribution  on  the  panel. 


Geometry 

The  first  step  In  panel  development  Is  the  definition  of  panel 
geometry.  Initially  assume  the  panel  is  a  trapezoid  lying  in  the  xy- 
plane  and  having  two  edges  parallel  to  the  x-axls .  It  is  then  sub¬ 
divided  into  two  triangles  having  a  common  side  that  joins  the  upper 
left  hand  corner  to  the  lover  right  hand  corner.  The  panel  is  oriented 
so  that  the  root  and  tip  chords  lie  parallel  to  the  free  stream  flow 
direction  at  Cf  *  0.  Figure  1  depicts  the  panel  geometry,  corner 
point  numbering  scheme,  and  coordinate  system. 


Singularity  Strength  Distribution 


The  general  form  of  the  singularity  strength  distribution  on  the 
panel  will  be 


Cl)(*.  y)  •  6(*.  y)T+  T(x,  y)7 


(2.1) 


where  and  "Y  are  continuous  functions  of  x  and  y.  For  the  purpose 
of  this  study,  (5  and  'Y  are  assumed  to  have  the  "bilinear"  forms 


<5(x.  y)  -  F  +  Dx  +  Ey 


(2.2) 


IfXx,  y)  -  A  +  Bx  +  Cy 


(2.3) 
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where  the  coefficients  A,  B,  C,  D,  F,  D  and  E  are  yet  to  be  determined 
constants.  These  constants  will  be  expressed  as  functions  of  the  panel 
geometry  and  the  unknown  singularity  strengths  at  the  corner  points. 

Consequence  of  the  Helmholtz  Condltldh 

The  vorticity  distribution  6l)(x.  j)  is  required  to  satisfy  the 
Helmholtz  condition  that  vorticity  must  be  preserved  in  the  fluid.  Thus* 

v<  +  Tj)  -  °  <2*4> 

on  the  panel.  If  §  and  y  have  the  toTUB  (2*2>  “d  (2*3>»  reBPec_ 
tively,  then 

0(F  +  Dx  +  Ey)/  3x  +  3(A  +  Bx  +  Cy)/  Qy  -  0  (2.5) 

which  implies 

D  -  -  C  (2.6) 

such  that 

5<x.  y)  -  F  -  Cx  +  Ey  (2.7) 

The  formulation  could  be  continued  in  terms  of  the  bilinear 
coefficients.  However,  it  is  conventional  to  express  them  In  terms  of 
panel  corner  point  vortidties. 

Bilinear  Coefficients  in  Terms  of  Corner  Vorticity 

In  the  ensuing  discussion,  the  subscript  L  refers  to  the  leading 
triangle  and  the  subscript  T  refers  to  the  trailing  triangle. 
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The  7  component  of  vorticity  (eq  (2.3))  Is  assigned  the  unknown 
▼•lues  yv  O3.  %  the  corners  (Xj,  y^,  (x3,  y2),  and  (x4,  y2) 
of  the  lending  triangle.  This  leads  to  the  following  system  of  three 

and  C^s 


equations  in  A^, 


'^'*i  +  *iIi  +  Vi  <2-e> 

X  -  *1 +  *i*j  +  V2  «•»> 

TJ  *  *L  +  *1%  +  V2  (2-10) 


The  system  has  the  solution: 


y2  ^V(y2  “  yl^  +  t(ylx4  "  xiy2^ /  (y2  ”  yi^ 

^/(x3  -  x^)  +  t (*xy2  -  y1x3)/(y2  -  y2)] 

1J/(x3  -  x4)  (2.11) 

»L  “  <  y3  -  T£>/<*3  -  *«>  <2*12> 

Cj,  -  -  TJ/(y2  "  yl*  +  I(xl  “  V/(y2  "  yl)J 

(x3  -  x4)  +  l(x3  -  x1)/(y2  -  y3)]  74/ 

(*3  -  x4)  (2.13) 


be  used  to  solve  for 
and  in  the  5  component  eq  (2.7).  Assigning  the  values  5^  and  5 3 
at  the  corners  (x^,  y^)  and  (x3,  y2)  leads  to 

61  +  Vi  •  Vi  +  h  <2-u) 


Since  -  - 


C^,  only  two  corner  conditions  may 
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+  V3  ‘  V2  +  rL 

and 

*J_  *  (  8x  -  63)/(7j  -  7j)  +  Klj  -  *3)/ 

(7j  -  Tj)!^ 

"  (yx  "  7 2  5^/(7j  “  y2)  ^  t ^lx3  ” 

Xly2)/(yl  "  y2>1CL 

In  a  similar  way,  the  trailing  triangle  coefficient  equations  are 
obtained  as  follows: 

Ay  ■  7l  ~  72^  ^  ^y2*l  ~  *4yl^^yl  "  ^2} ^ 

7 f2/<*2  “  xl^  +  ^*4yl  “  y2x251  7j/ 

(*2  -  *!> 

By  -  (  Y2  ~  7j)  /  (X2  "  *1^ 

Cy  -  -  -  72>  +  I(*4  -  *£Kj-y  -  y2)l  y  1! 

(*2  -  *j)  ♦  l(*2  “  *4)/(7i  “  y2)l  Tj/ 

<x2  "  xl) 

*T  "  <  52  "  4>/<yl  '  y2>  +  I(x2  "  xk>/ 

(yl  “  72)3ct 


(2.15) 

(2.16) 

(2.17) 


(2.18) 

(2.19) 

(2.20) 

(2.21) 
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FT  -  (7i  §4  ”  y2  52>/<yi  "  V  *  1^1*4  " 

~  y2)lcT  (2.22) 

where  and  have  been  expressed  in  terms  of  (as  functions  of )  (j) 
components  at  corners  2  and  4. 

Mathematical  Continuity 

The  bilinear  vorticity  distribution  is  continuous  on  each  triangular 
region.  In  addition,  the  'Y  component  has  been  made  to  be  continuous 
throughout  the  planform  by  the  representation  In  terms  of  corner  values. 
This  can  be  demonstrated  as  follows.  Let  be  the  boundary  shared  by 
any  two  adjacent  triangles.  Then  J-1  Is  a  straight  line  segment  and  can 
be  deecribed  by  a  linear  expression  (i.e.,  y  In  terms  of  x  or  x  in  terms 
of  y).  The  'Y  distribution  on  each  of  the  adjacent  triangles  will 
degenerate  to  a  linear  function  of  a  single  variable  upon  substitution 
of  this  expression.  Both  functions  assume  the  same  "Y  values  at  the 
endpoints  of  J”*  .  Since  only  two  points  are  needed  to  determine  a 
straight  line  or  a  linear  form,  we  have  'Y  watching  identically  on 

r. 

The  §  component  has  breaks  in  continuity  throughout  the  planform. 
This  is  a  consequence  of  applying  the  Helmholtz  condition  (eq  (2.4)) 
which  eliminated  the  constant  D  and  expressing  the  remaining  two 
unknowns  in  terms  of  two  $  corner  values,  out  of  a  possible  three. 

The  (5  component  is  continuous  on  panel  leading  and  trailing  edges 
since  it  is  on  these  edges  that  conmon  (5  values  are  assumed  at  the 
endpointe.  The  $  component  is  discontinuous  on  panel  diagonals. 
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Application  of  the  Blot-Savart  Law 

Let  CO  (f)  be  a  vorticity  distribution  defined  on  a  finite  region 
R  in  the  xy  plane.  Let  $  be  a  fixed  point  (control  point)  In  the  plane. 
The  velocity  Induced  at  s  due  to  the  distribution  on  R  Is  given  by  the 
Biot-Savart  Lav  (Ref  5:526-528): 


47Tv(*) 


f 


~>3. 


[  CJ(*)  x  (s  -  *)]/|S  -  x|dR 


(2.23) 


Suppose  the  control  point  s'  la  located  at  the  origin  of  the 
coordinate  system.  Note  that  this  can  be  done  by  performing  a  simple 
translation  of  the  plane.  Then, 


s  -  (0,  0,  0)  (2.24) 

t  -  *  -  (0,  0,  0)  -  (x,  y,  0)  -  (-x,  -y,  0)  (2.25) 


Assuming  the  distribution  has  the  form  (2.1), 

G)X  (f  -  f)  -  1x7  ~  y$  Jk  (2.26) 

vhere  C  Is  the  unit  vector  normal  to  the  xy-plane.  Substituting  the 
expressions  (2.25)  and  2.26)  Into  eq  (2.23)  yields: 


47TV  (0,  0,  of 


-  7  d  >/<*2  +  y2>3/2« 


(2.27) 


vhere  w  is  the  normal  velocity  component  Induced  at  the  origin. 
Substituting  the  expressions  for  7  (2.3)  and  $  (2.7)  Into  eq  (2.27) 

yields: 
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n.-JJ 


(Ax  +  Bx2  +  2Cxy  -  EyZ  -  Py)/ 


(,2  ♦  jh3,1a 


where  R  is  taken  as  the  region  defined  by  a  trapezoidal  panel. 


Let  and  R^  be  the  subregions  of  R  which  correspond  to 
Ing  and  trailing  triangles.  The  coefficients  A,  B,  C,  F,  and 
constant  on  each  subregion  and  eq  (2.28)  is  rewritten  as: 

47T—  Vi  +  vl  *  V2  +  V*  +  2  Vs  + 
*£  '  *£  -  V*  *  Vs  -  Vs 

where 

T  -  II  (x/(*2  +  y2)3/2)dR 
*L,  T 

I2>  T  “  if  (*2/<*2  +  *2>3/2>dR 

V  T 

i\*  T  •  //  (x^/<5t2  +  y2>3/2>dR 

*L.  T 

ij*  T  -  //  +  y2)3/2>dR 

\  T 

Ij*  T  •  J I  (y2/(*2  +  y2)3/2)dR 
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(2.28) 

the  lead- 
E  remain 


(2.29) 

(2.30) 

(2.31) 

(2.32) 

(2.33) 

(2.36) 
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Evaluation  of  these  integrals  may  be  found  in  Appendix  A.  Substitution 
of  the  expressions  for  A^  through  and  collecting  coefficients  of  the 
unknowns  (  Sy  52*  53*  64.  0^.  "Yy  and  O4)  yields,  after 

considerable  algebraic  manipulation, 

4  7T  W  -  C(y  jY  i1}/^  -  y2)l  6X  + 

r(y2x4  ‘  -  y2)]  52  + 

I(I5“ylI4)/(yl"y2)1  ^3  + 
t(15  "  ylX4)/(yl  "  y2)]  ^4  + 

[((*47!  “  y 2*2^1  ”  ^x2  "  *i)y2yi  “ 

(yl  -  72>x2  +  ^*2  ”  *1^*1  + 

(x2  "  X4>V/((X2  "  xl)(yl  "  y2))]  Ti  + 

[((y2xi  “  +  (Ji  ~  72)X2  + 

(*4  “  *1)KT)/((y1  -  72)(x2  “  *!»]  *Y2  + 
t((ylx4  “  x172)I1  +  ^2  “  yl)X2  + 

<xx  “  x4)\)/«y2  “  TjKxj  -  X4))]  %  + 

t((Vz  -  V3»l  -  <x3  -  Wl  - 

(y2  -  7j)l\  +  <x3  -  Xj)\  +  (x3  -  x4)Kt>/ 
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((y2  -  -  *4)>1  T4  (2.35) 

where 

Kj,  “  21^  -  [(x1  -  *3)/^  -  y2))I^  - 

1(72*3  “  *172^^1  “  ^2^T4  (2.36) 

*1  -  21*-  IC*2  -  *4)/^  -  72)llJ  - 

f  <73*4  -  *272)/(72  -  72)]lJ  (2.37) 


Summary 

Expression  (2.35)  Is  the  normal  velocity  Induced  at  the  origin  of 
the  xy  plane  by  a  trapezoidal  vortlcity  panel.  This  velocity  Is  due  to  a 
bilinear  vortlcity  distribution  which  satisfies  the  Helmholtz  condition 
eq  (2.4). 
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III.  Panel  Assembly 

This  section  presents  the  panel  assembly  procedures  needed  to 
model  a  wing.  The  goal  Is  to  develop  the  methodology  required  to 
predict  the  pressure  distribution  and  associated  forces  and  moments  on 
a  wing. 

Panel  Numbering 

Figure  2  Illustrates  a  paneling  arrangement  and  associated  number¬ 
ing  scheme  for  a  16  panel  vlng.  The  panels  are  numbered  consecutively 
In  the  chordwise  direction  starting  with  the  inboard  leading  edge 
panels  and  terminating  tr*tl.  the  outboard  trailing  edge  panels. 

Number  of  Unknowns.  Boundary  Conditions  and  numbering 

Let  M  be  the  number  of  chordwise  panels  and  N  be  the  number  of 
spanwlse  panels.  These  are  defined  using  M  +  1  chordwise  cuts  and 
N  +  1  spanwlse  cuts.  Each  Intersection  determines  a  panel  corner 
point.  Since  there  ere  two  unknown  components  at  each  corner  point, 
the  total  number  of  unknowns  is  given  by: 

2  (M  +  1)  (N  +  1)  (3.1) 

Boundary  Conditions 

TVo  boundary  conditions  are  imposed  on  the  wing  panel  system. 

These  reduce  the  number  of  unknowns  and  improve  the  physical  modeling 
of  the  flow  field. 
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The  z  axis  1b  normal  to  the 
wing  plan form. 


X 


Figure  2.  A  Paneling  Arrangement  and  Associated 
Numbering  Scheme  for  a  16  Panel  Wing 
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The  Kutta  condition  (Ref  5:390-399) 

7(x,  y)  -  0  (3.2) 

is  imposed  at  all  comer  points  on  the  ving  trailing  and  tip  edges. 

Since  the  ~Y  component  of  vorticity  eq  (2.2)  is  both  linear  and 
continuous  on  these  edges,  the  Kutta  condition  is  satisfied  Identically. 
The  Kutta  condition  reduces  the  total  number  of  unknowns  by: 

M  +  N  +  1  (3.3) 

(NOTE:  The  point  defining  the  Intersection  of  the  trailing  and  tip 
edges  is  common  to  both.) 

A  historically  acceptable  boundary  condition  (an  outgrowth  of 
Prandtl's  lifting-line  theory  (Ref  5:535-567))  is  for  the  vorticity 
vector  to  lie  tangent  to  the  wing  leading  edge.  This  boundary  condi¬ 
tion  initially  orients  the  vorticity  vector  so  that  a  positive  circu¬ 
lation  is  produced.  The  boundary  condition  is  imposed  at  all  leading 
edge  comer  points  and  can  be  written: 

7/<5  -  A  (3.4) 

or 

S  -  7/A  (3.5) 

where  A  is  the  leading  edge  slope  at  the  comer  point.  It  reduces 
the  total  number  of  unknowns  by: 

N  +  1  (3.6) 
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The  unknown  corner  n  a  are  denoted 
by  **-  ”  and  the  7  s  by  "■ 


Note:  7-0  ♦  0 

at  the  'LE  tip  ^ 
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and  the  total  number  of  unknowns  becomes: 

2MN  +  M  (3.7) 

after  imposing  the  two  boundary  condition  equations  (3.2)  and  (3.5). 

Since  the  wing  is  symmetric  about  the  z-axls,  it  may  seem  logical 
to  Impose  a  boundary  condition  at  the  wing  root  chord.  However,  setting 

5“  0  (3.8) 

at  the  centerline  is  redundant  for  rectangular  wings  and  leads  to  an 
ill-conditioned  system  once  planform  symmetry  is  considered. 

Unknown  Numbering  Scheme 

Figure  3  illustrates  the  unknown  numbering  scheme  for  a  9  panel 
wing  with  applied  boundary  conditions.  The  paneling  arrangement  of 
Figure  3  is  chosen  because  it  represents  the  smallest  number  of  panels 
needed  to  illustrate  interior  panels  and  panels  having  boundary 


conditions.  Let  1  be  the  panel  number  of  an  interior  panel.  Then 

the  following  numbers  (in  terms  of  1,  M,  and  N)  are  assigned  to  the 

eight  unknown  corner  vortldty  components: 

(6X) 

i  - 1 

(3.9) 

(fi2) 

i 

(3.10) 

(53) 

i  +  M  -  1 

(3.11) 

<54> 

i  +  M 

(3.12) 
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<Ti> 

MN  +  M  +  i 

(3.13) 

<r2> 

MN  +  M  +  i  +  1 

(3.14) 

<V 

MN  +  2M  +  i 

(3.15) 

<X> 

MN  +  2M  +  i  +  1 

(3.16) 

Solving  for  the  Corner  Vorticitles 

The  total  number  of  unknowns  (after  the  boundary  conditions  are 
applied)  is  given  by  eq  (3.7)  which  also  specifies  the  number  of  condi¬ 
tions  needed  to  solve  for  the  corner  vorticities.  Two  types  of  condi¬ 
tions  will  be  used;  control  point  conditions  and  edge  or  point 
continuity  conditions. 

Control  Point  Equations 

Control  point  equations  are  obtained  using  eq  (2.35).  The 
velocity  component  is  computed  for  one  point  (control  point)  on  each 
panel  comprising  the  wing. 

Let  (x^,  y^)  be  the  control  point  on  panel  i.  To  obtain  the 
contribution  to  w^  due  to  panel  j,  express  the  coordinates  of  panel  j 
in  a  coordinate  system  with  (x^,  y^)  at  the  origin.  This  is  done  by 
performing  a  translation  in  the  z  ■  o  plane: 


Equation  (2.35)  is  then  applied  with  appropriate  boundary  conditions. 
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Planform  symmetry  Is  Included  by  reflecting  either  the  panel  or  con¬ 
trol  point  about  the  x-axls  and  applying  eq  (2.35).  Reflecting  the 
control  point  Is  less  complicated  from  a  programming  viewpoint. 


The  above  process  is  repeated  for  each  panel  on  the  vlng.  After 


ell  the  contributions  to  have  been  calculated,  it  can  be  written  as: 


2MN  +  M 

I 

i  -  l 


k±A 


(3.18) 


where  is  the  column  vector  of  unknown  corner  vortlcltles  and  the 
coefficients  A^  are  functions  of  panel  geometry.  The  are  numbered 
using  the  system  given  by  eqs  (3.9)  through  (3.16).  One  control  point 
equation  (3.18)  Is  obtained  for  each  panel  on  the  wing  and  together 
they  comprise  MN  conditions. 


Edge  Continuity  Conditions 

The  §  component  of  vorticity  eq  (2.7)  Is  discontinuous  across 
the  panel  diagonal  (See  discussion  In  Section  II.).  This  can  be 
partially  remedied  by  specifying  a  point  continuity  condition  at  the 
panel  lower  right  hand  corner.  This  condition  Is: 

Sl («4.  y1i  -  4  o.i9) 

which  becomes 


"63  +  ^x3  x4^y2  ”  71^  ^1  “ 

[(xj^  -  x4)/(y2  “  7i)3  %  "  I(*3  "  *j>/ 

(y2  '  7i>J  X  "  0  0.20) 


20 


AFIT/GAE/AE/ 80J-1 


after  substituting  of  (x^,  y2)  end  the  expressions  for  A^  (2.11). 

(2.12).  end  (2.13)  Into  eq  (2.7).  One  condition  eq  (3.20)  is 
formulated  for  each  panel  on  the  wing  for  a  total  of  MN  conditions. 

Note  that  £)  is  still  not  continuous  across  the  panel  diagonal  due 
to  the  remaining  discontinuity  at  the  upper  left  hand  corner.  Also. 

1)  is  not  continuous  across  panel  side  edges. 

The  edge  continuity  conditions  and  control  point  equations  total 
2MN  conditions.  M  additional  conditions  can  be  obtained  by  specifying 

an  edge  continuity  condition  for  §  at  the  upper  left  hand  corner  of 

\ 

each  panel  along  the  centerline.  This  condition  is: 

&<*!»  yx)  -  (3.21) 

which  becomes 

6-2  -  <5i  -  iu2  -  xi>/(y!  -  + 

[(*4  ~  “  7^y2  +  ( (*2  ” 

(yx  -  y2)  Ti  •  0  O.22) 

after  substitution  of  (x^,  y^)  and  the  expressions  for  A^  (2.18). 

(2.19).  and  (2.20)  into  eq  (2.7).  This  choice  is  based  on  trial  and 
error,  the  additional  (5  continuity  on  the  centerline  having  the  effect 
of  minimizing  vorticity  oscillations. 

Compressibility 

Compressibility  effects  are  accounted  for  by  using  the  Prandtl- 
Glauert  transformation  (Ref  2:124-127): 
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x  -  x/  -yjl  -  M2  (3.23) 

The  transformation  la  applied  to  all  x  coordinates  which  are  used  in 
either  the  control  point  equations  or  edge  continuity  conditions. 


Matrix  Formulation 


The  control  point  equations  and  edge  continuity  conditions  are 
2MN  +  M  equations  in  the  unknowns,  Q  .  This  system  has  the  matrix 
formulation: 


-Gl- 

’  W1  ' 

[Aj  * 

0 

-  ' 

''MN 

2MN  +  M 

o 

(3.24) 


(2MN  +  M)x(2MN  +  M) 


The  first  MN  rows  of  are  the  coefficients  fcr  the  control  point 
equations  and  are  all  nonzero.  The  last  MN  +  M  rows  of  correspond 
to  the  homogeneous  edge  continuity  conditions  and  have  no  more  than 
five  nonzero  entries  per  row. 


Solution 

The  linearized  form  of  the  flow  tangency  boundary  condition  is 
(Ref  5:495): 


w/voo  dc/d*  ~Ct  (3.25) 

where  Qt  the  wing  angle  of  attack,  dc/dx  is  the  local  camber  slope 
and  V^is  the  free  stream  velocity.  This  expression  is  substituted 
for  each  w^  in  eq  3.24  where  dc^/dx  is  the  panel  slope  at  control 
point  i.  In  matrix  notation, 
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MM 


Wdci/dx  ~a) 

o 


(3.26) 


or 


r.r\ 

dct/dx  -Qf 

[  uj  • 

o 

(3.27) 


Forces  and  Moments 

Once  the  ^oo  are  °^ta^ne^  for  a  given  Qf  **6  camber  slope 
distribution,  eqs  (2.2)  and  (2.7)  can  be  used  to  calculate  the 
vorticlty  strength  at  any  point  on  the  planform.  The  surface  perturba¬ 
tion  velocities  In  terms  of  local  vorticlty  are  (Ref  5:508): 


“^co  • 

±  T/2 

(3.28) 

+  6/2 

(3.29) 

The  upper  sign  corresponds  to  the  upper  wing  surface  and  visa-versa. 

Pressure  coefficients  are  obtained  from  the  perturbation  velocities 
by  either  using  the  exact  isentroplc  expression  (Ref  3:167): 

Cp  -  2[(1  +  (7-  1)M^/2(1  -  (0^+  u)2  + 

v2  +  w2)  /V^)  y,{y~  15  -  11/(7^)  (3.30) 

or  the  second  order  approximation  (Ref  3:167): 

Cp  -  -  [2u/V00+  (1  -  M^>)u2/V00+  (v2  +  v2)/V^]  (3.31) 

which  Is  adequate  for  two-dimensional  and  planar  flows.  These  coeffi¬ 
cients  are  Integrated  along  chord  lines  to  obtain  local  lift  and  moment 
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coefficients. 


The  appropriate  expressions  are 


«t 


Si 


C  )dx 
P.. 


C  )x  dx 
P„ 


(3.32) 


(3.33) 


where  the  subscripts  1  and  u  refer  to  the  lover  and  upper  wing 
surf aces . 

This  concludes  the  theoretical  section  of  this  report.  The  next 
step  Is  the  Implementation  of  this  panel  to  predict  airloads  on  vlnge. 
This  is  done  by  the  use  of  a  computer  code  "WING”  which  is  presented  in 
Section  IV. 
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IV.  Computer  Code 


General  Description 

A  FORTRAN  code  "WING”  has  been  developed  to  analyze  planform  flow 
using  the  methodology  discussed  In  Sections  II  and  III.  WING  la  a 
pilot  code  and  should  not  be  treated  as  a  fully  checked-out  production 
program  until  the  problems  In  the  wing  tip  region  are  resolved 
(Section  V).  Many  of  the  programming  techniques  used  In  WING  have 
been  previously  developed  by  the  author  and  can  be  found  In  Reference  8. 
Care  has  been  taken  to  Insure  correspondence  between  FORTRAN  variable 
names  and  the  symbol  usage  in  Sections  II  and  III.  The  listing 
(Appendix  C)  contains  comment  cards  that  describe  program  function  and 
logic  In  detail.  Below  Is  a  brief  description  of  each  of  the  sub¬ 
routines  In  WING. 


WING  (Main) 

WING  Is  the  executive  control  routine.  All  geometric  data  Is 
read  by  WING.  WING  Initializes  panel  parameters  and  calls  subroutines 
MESH,  AERO,  INVRT,  and  LOADS  in  that  order. 

MESH 

MESH  generates  the  x  and  y  coordinates  for  the  panel  corner  points 
and  control  points.  The  mesh  Is  generated  from  the  Information  given  on 
the  first  four  data  cards.  MESH  Is  a  FORTRAN  version  of  the  mesh 
generator  discussed  In  Reference  8. 

25 
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AERO 

AERO  formulates  the  control  point  equations  and  edge  continuity 
conditions.  Subroutine  INT  is  called  by  AERO. 

INT 

INT  evaluates  the  panel  Integrals  using  equations  (A. 31)  and 
(A. 32).  Subroutine  STRIP  Is  called  by  INT. 

STRIP 

STRIP  evaluates  the  five  STRIP  function  equations  (A. 24)  through 
(A. 28)  given  two  points  in  the  plane. 

INVRT 


INVRT  Inverts  the  coefficient  array  for  the  system  of  equations 
formulated  by  AERO.  The  inversion  is  performed  using  Gaussian 
elimination.  INVRT  la  essentially  the  same  Inversion  subroutine  found 
in  the  7ASTL0DS  lifting  surface  program  (Ref  6:73-133). 

LOADS 

LOADS  calculates  the  planform  pressure  distributions  and  aero¬ 
dynamic  coefficients  for  the  loading  cases  specified  by  cards  6  and  7. 
After  the  last  loading  case  is  examined,  LOADS  will  terminate  program 
execution. 

Input  Description 

This  section  provides  a  card  by  card  description  of  Input  data 
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•long  with  some  helpful  "dos  and  don'ts"  of  program  operation.  All 
input  la  unformatted  and  should  be  separated  by  commas.  Integer  data 
cannot  have  a  decimal  point.  WING  la  a  nondlmenslonal  code  which  will 
accept  data  In  any  consistent  system  of  units.  Presently,  WING  can 
analyze  plan forms  having  sixty  panels  or  less.  This  can  be  Increased 
by  changing  the  dimensions  of  arrays  X,  Y,  XC,  YC,  E,  A,  SG,  CBR,  and 
SUM  in  common  blocks  A  and  C.  The  reader  is  referred  to  Appendix  B 
which  contains  a  sample  problem. 


Card  1  (Span  Data) 


Variables  (In  Order) 
SSPN 

SS 


S(l),  ...,  S (NS) 


Description 

The  length  of  the  wing  semi-span. 

The  number  of  stations  needed  to 
define  the  spanwlse  panel  boundaries. 
The  wing  root  chord  is  station  1  and 
the  last  station  1b  the  wing  tip.  NS 
is  an  Integer  and  must  be  less  than  or 
equal  to  12. 

Span  stations  as  a  fraction  of  the 
semi -span.  0.  and  1.  will  always  be 
the  first  and  last  entries.  Entries 
must  be  In  ascending  order  left  to 
right. 
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Card  2  (Chord  Data) 

HC  The  number  of  stations  needed  to 

define  the  chordvise  panel  boundaries 
The  leading  edge  is  station  1  and  the 
last  station  is  the  trailing  edge.  NC 
is  an  Integer  and  must  be  less  than  or 
equal  to  10. 

C(l),  . ...  C(NC)  Chord  stations  as  a  fraction  of  the 

wing  local  chord.  0.  and  1.  will 
always  be  the  first  and  last  entries. 
Entries  must  be  in  ascending  order 
left  to  right. 

Card  3  (Break  Point  Data) 

Break  points  are  the  z  coordinates  of  leading  and  trailing  edge 
for  those  chords  that  define  a  sweep  change.  The  wing  centerline 
defines  the  positive  z-axls  with  origin  at  the  leading  edge. 

NB  The  number  of  break  point  sets 

needed  to  outline  the  plan form 
geometry.  Two  sets  will  be  needed 
to  define  a  four-point  wing.  NB  is 
an  integer  and  must  be  less  than  or 
equal  to  10. 
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NL(1) ,  B(l,  1),  B(l,  2),  ....  NL  is  the  number  (Not  the  value!)  of 
NL(NB) ,  B(NB,  1),  B(NB,  2)  the  span  station  vhere  the  sweep  change 

occurs.  NL  is  an  Integer.  B(I»  1)  and 
B(It  2)  are  the  leading  and  trailing 
edge  x  coordinates  of  the  chord  line 
at  span  station  NL(I).  Break  point 
sets  are  entered  from  Inboard  to  out¬ 
board.  The  first  set  (NL  *  1)  will 
always  he  the  x  coordinates  of  the 
root  chord  leading  and  trailing  edge. 
The  last  set  (NL  ■  NS)  will  always  be 
the  x  coordinates  of  the  tip  chord 
leading  and  trailing  edge. 

Card  4  (Mach  and  Control  Point  Data) 


CY,  CX  Control  point  location  in  terms  of 

local  panel  span  and  chord.  CY  is 
the  fraction  of  the  panel  span  and  CX 
is  the  fraction  of  the  panel  chord. 

A  recommended  control  point  choice  16 
CY  ■  .15  and  CX  •  .75  which  is  based  on 
extensive  program  testing. 

MACH  Mach  number  (not  an  integer) .  WING 

accepts  subsonic  Mach  numbers  only. 
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Card 

NA 


NA  cards  are  required 

ALPHA 

NCCG 


NPF 


(Number  of  Loading  Cases) 

The  number  of  loading  cases  to  be 
examined  -  1. e . ,  changes  in  angle  of 
attack  or  camber  slope  distribution. 
NA  Is  an  Integer  that  has  no  upper 
bound. 

Card  6  (Load  Cases) 


The  angle  of  attack  In  degrees.  Plus 
Is  nose  up. 

The  camber  change  parameter  (integer) . 
Enter  0  to  read  a  new  camber  slope 
distribution  or  1  to  retain  the  pre¬ 
vious  distribution.  WING  Initially 
sets  the  camber  slope  array  equal  to 
0. 

The  pressure  option  parameter 
(Integer) .  Enter  0  to  use  the  exact 
lsentropic  expression  eq  (3.30)  or 
enter  1  to  use  the  2nd  order  approxi¬ 
mation  eq  (3.31). 
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CY,  CX 


This 

CBR(l)  .. 


Panel  location  (fraction  of  panel 
spas,  fraction  of  panel  chord)  where 
the  pressures  are  to  be  calculated. 
Pressures  may  be  computed  at  points 
other  than  the  control  points. 

Card  7  (Camber  Slope  Distribution) 
card  is  used  only  for  a  NCCG  value  of  0. 

,  CBR(NP)  The  value  of  the  local  panel  slope  In 

degrees  for  each  panel  on  the  plan- 
form.  NP  is  the  total  number  of 
panels.  Entries  must  be  made  in  the 
order  corresponding  to  the  panel 
numbering  scheme  (Figure  2).  The 
sign  rule  for  camber  follows  the 
standard  convention. 
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V.  Results 

Program  WING  was  exercised  for  a  variety  of  four  point  wings 
having  various  aspect  ratios,  taper  ratios,  and  sweep  angles.  This 
section  presents  results  for  two  of  these  wings. 

Two  general  observations  are  made  first.  One,  control  point 
location  is  the  major  factor  controllng  bounded  numerical  oscillations 
of  the  vortldty  vector  as  it  is  in  many  current  paneling  routines 
(ex.  Refs  6  and  8) .  Oscillations  are  very  common  if  the  control 
point  is  located  anywhere  on  the  leading  triangle.  Fewer  oscillations 
occur  if  the  control  point  is  located  on  the  trailing  triangle  with 
.1  *  CY  *  .5  and  .A  •*  CX  -  .9.  Secondly,  the  program  shows  the 

best  results  when  uniform  spanwlse  paneling  is  used.  Non-uniform 
span wise  paneling  tends  to  cause  oscillations  in  the  vortlclty  vector. 
However,  non-uniform  chordwlse  paneling  seems  to  have  little  effect  on 
solution  stability.  The  best  total  match  (with  other  known  solu¬ 
tions)  occurs  at  approximately  CY  ■  .15  and  CX  •  .75. 

Rectangular  Wing 

The  first  case  examined  is  a  rectangular  wing;  AR  ■  8,  Ct  ■  5°  and 
Mqq  ■  .1.  The  wing  is  modeled  using  12  uniformly  spaced  span  stations 
and  6  non-unlformly  spaced  chord  stations  (0.,  .1,  .3,  .6,  .8,  and  1.) 
which  define  55  panels.  Figure  A  shows  the  distribution  predicted  by 
WING  and  Anderson  (Ref  1:9-16).  The  lift  coefficient  experiences  a 
spiking  phenomena  near  the  wing  tip.  This  phenomena  always  happens  in 


AFIT/GAE/AE/80J-1 


the  region  defined  by  the  last  three  outboard  span  stations.  Also,  the 
distribution  falls  off  faster  than  it  should  before  reaching  the 
region  where  spiking  occurs.  There  is  very  little  difference  between 
total  lift  coefficients  indicating  the  areas  under  the  curves  are 
approximately  equal.  This  shows  the  solution  is  possibly  trying  to 
compensate  for  the  spike  by  underpredicting  lift  in  the  Inboard  region. 

Figure  5  shows  the  spanwlse  distribution  of  center  of  pressure. 

The  Xc  shifts  are  aft  in  the  region  of  spiking  before  traveling  for¬ 
ward.  The  ACp  versus  chord  station  curves  of  Figure  6  exhibit 
expected  behavior  for  the  inboard  stations  (0  *  Y/S  *  .585).  At 
Y/S  -  .9,  the  curve  has  "fattened  up"  considerably  which  drives  the 

X  backwards  and  creates  the  C.  spike.  The  curve  is  subsiding  again  at 
cp  L 

Y/S  -  .95  since  the  Kutta  (no  net  load)  boundary  condition  is  imposed  at 
the  tip. 

Figure  7  is  the  7  strength  at  the  root  chord.  The  calculated 
solution  compares  favorably  with  the  exact  2-D  flat  plate  solution 
(Ref  5:515)s 

7(X)  -  2a  [(C  -  x)/(cx  -  x2)]  (5.1) 

where  c  is  the  chord  length  and  Q f  is  measured  in  radius. 

Figures  8  and  9  show  7  ®ad  (5  strength  distributions  at 
selected  span  stations.  The  £  distribution  grows  in  magnitude  rela¬ 
tive  to  7  *s  ve  approach  either  the  tip  or  the  trailing  edge.  This 
allows  the  vorticlty  vector  to  turn  and  satisfy  the  Kutta  condition  as 
shown  in  Figure  10.  The  (5  component  is  0  at  the  leading  edge  which 
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ACP 


Figure  6.  Versus  X/c  at  Selected  Span  Stations 

Rectangular  Wing,  AR  -  8,  Ct  ■  5° 
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is  the  condition  of  infinite  slope  (eq  (3.5)).  Also,  the  ($  component 
is  small  near  the  wing  root  chord  which  illustrates  the  effect  of  plan- 
form  symmetry.  Both  the  §  and  'Y  components  exhibit  unusual 
behavior  in  the  tip  region.  Possible  remedies  for  this  problem  are 
discussed  in  Section  VI. 

Swept  Wing 

The  second  case  examined  is  a  swept  un tapered  wing;  AR  «  4.5, 

A*  40°,  Ct  -  5°,  and  “  .1.  The  wing  modeling  is  the  same  as 

the  rectangular  wing.  Figure  11  shows  the  distribution  predicted  by 
WING  and  by  wind  tunnel  tests  (Ref  4:92).  Again,  a  spiking  phenomena 
occurs.  The  predicted  curve  16  showing  the  proper  curvature  in  the 
region  Inboard  of  the  spike.  Also,  the  total  lift  coefficients  again 
show  close  agreement.  Figure  12  shows  the  x  center  of  pressure  versus 
span  station  which  again  shifts  aft  in  the  region  of  spiking. 
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is  Span  Station  for  a  Sveptback 
id  King;  AR  -  4.5,  A  -  40°,  Ct  -  5° 
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71.  Conclusions  and  Recommendations 


Conclusions 

A  paneling  technique  vhich  allows  the  vorticlty  vector  to  change 
direction  has  been  demonstrated.  The  proper  choice  of  control  point 
(see  Section  V)  will  guarantee  a  solution  free  of  numerical  oscillations 
in  the  vorticlty  vector.  This  method  has  been  Implemented  on  a  computer 
and  Introduces  no  new  complexities  to  an  experienced  programmer.  Exist¬ 
ing  mesh  generators  and  other  aerodynamic  modules  were  incorporated  into 
this  technique  (Ref  8) .  Computer  run  times  are  of  the  same  order  as 
programs  incorporating  "fixed  direction  vorticlty"  panels  and  no 
problems  involving  extensive  run  times"  vere  encountered. 

The  method  gives  good  aerodynamic  results  near  the  centerline  of 
the  wing.  However,  the  method  will  underpredict  life  as  we  move 
outboard.  A  gross  overprediction  of  lift  occurs  In  the  region  defined 
by  the  three  most  outboard  span  stations.  Total  lift  coefficients  as 
predicted  by  this  technique  agree  well  with  existing  solutions. 

Recommendations 

The  following  ideas  are  suggested  for  the  Improvement  of  the 
method  and  hopefully  will  lead  to  the  elimination  or  minimization  of 
the  "spiking"  problem. 

a.  A  wake  model  should  be  incorporated  into  the  program. 

b.  The  flat  plate  panels  should  have  a  provision  for  leading  edge 
suction . 
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c.  Higher-order  panels  might  be  needed  at  the  leading  edge  and 
tip.  This  would  turn  the  vorticity  vector  faster  and  minimise  the 
spiking.  An  elliptic  vorticity  distribution  is  suggested  since  many 
classic  lift  distributions  are  elliptic  near  the  tip. 

d.  Change  or  modify  the  leading  edge  boundary  condition.  The 
nclasslcal''  tangency  condition  may  be  unappropriate  for  this  kind  of 
panel.  Reference  9  suggests  the  Kutta  condition  be  applied  at  the 
leading  edge. 

e.  Interchange  the  role  of  “Y  and  with  respect  to  panel 
boundary  continuity  conditions.  Perhaps  the  “Y  distribution  enjoys 
too  much  continuity  on  the  planform.  This  could  be  creating  problems 
by  forcing  the  “Y  distribution  to  undergo  unusual  "warping"  in  order 
to  satisfy  the  planform  boundary  conditions. 
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and 


f3(x.  y)  -  xy/<*2  +  y2)3/2 
*4(x,  y)  -  y/(x2  +  y2)3^2 


f5(x,  y)  -  y2/(x2  +  y2)3^2 


(A.  7) 
(A.  8) 


(A.  9) 


To  Illustrate  the  methodology  used  in  evaluating  the  Integrals  (A.  4), 
consider 

rT  i 

*»  <)  2^/2 

+  yV'  ]dxdy  (A. 10) 

h)/M 


u*>  f 1  f 

L-*OCT  y  (y  -  1 


By  Fierce's  integral  tables  (Ref  7), 


F5  -  11a 


L*00  y 


|Tl  trV(7 


]dy  -  (A.  11) 


(y  -  b)/M 


11m 


L*CO  y 


f  1  L  dy/Vh2  +  y2  -  J  1  ((y  -  b)/M)dy/ 

Z  W  V 


Vc(y  -  W/M)2  +  y2 


(A.  12) 


Substituting 


*  -  (y  -  b)/K 


(A.  13) 


and 


dy  “  Mdz 


Into  the  rightmost  Integral  of  (A. 12)  yields: 


(A. 14) 


49 


AEIT/GAE/AE/ 80 J-l 


11m 

L*CO 


f 


(Mz  +  b)2 


(A.  15) 


The  Integrals  (A. 15)  can  be  evaluated  using  Ref  7.  One  obtains,  after 
some  algebraic  manipulation, 

P5  -  11m  L  ln[  (yx  +  VyJ  +  L2)/(yo  +  +  L2)]  + 

L+OO 

[b^/a  +  M2)3/2]ln[(\Cf+^2  +  ((xj  +  Hy^/ 

V1  +  *?))/(  Vxo  +  yo +  ((xo +  My0)/V1  +  m2))]  + 

[M/a  +  m^kV*2  +  yj  *  *V^i  +  yi>  ca.i6) 

recalling 

yo  -  Mxo  +  b  (A.  17) 

and 


yx  -  Hx1  +  b 
Notice  the  limit 


(A.  18) 


11m  L  ln[  (yx  +  ^f+T2>/(70  +  ^  +  L2)  ]  (A.  19) 

L-*CO 

is  not  a  function  of  x  or  x,  which  leads  to  the  following  observation. 

o  1 
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If  the  Integral  F^  le  evaluated  on  any  other  semi- infinite  atrip  bounded 
by  the  lines  y  ■  yQ  and  y  “  y^,  the  limit  (A. 19)  is  Invariant.  Evalu¬ 
ating  the  integral  F^  using  any  other  line  segment  connecting  y  ■  yQ  to 
y  -  y^  and  formulating  the  difference  betveen  this  result  and  (A. 16) 
leads  to  cancellation  of  the  limit  (A. 19). 

The  other  four  integrals,  Fj  through  F^,  can  be  evaluated  by  a 
similar  procedure  using  Ref  7.  Each  Integral  has  a  limit  term  given 
by: 


(Fl) 

lim  ln[(yo  +  Vy2  +  I2)/^  +  ^h\  +  l*2)) 
L-*CO 

(A.  20) 

(f2) 

lim[yiln (Vh2  +  j\  +  L)  -  yoln(Vh2  +  y2  + 

LvCO 

h)l 

(A.  21) 

(F3) 

lim[  +  y*  -  Vh2  +  y2l 

L-*CO 

(A. 22) 

<v 

lim  in[(L  +  Vl2  +  y2)/a  +  Vl2  +  yj)l 
L-*CO 

(A.  23) 

These  terms  cancel  upon  the  formulation  of  integral  differences. 
Retaining  the  finite  terms  from  each  of  the  Integral  evaluations,  ve 
can  define  "Strip  Functions"  for  the  functions  f^  and  the  points 

<v  yc>  (xi»  yi> by 

Sl[(xo*  yo}*  (xl*  yl)]  "  (M/Vl  +  ?)G  (A. 24) 
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S2I(V  V*  (X1‘  71)3  "  +  *0l  - 

+  F*  +  Xjl  +  (M/(l  +  M2))H  + 

[b/(l  +  M2)3/2]G  (A. 25) 


S3C(xo,  yo),  (Xj,  7l)J  -  (1^/(1  +  M2))H  + 

[bM/(l  +  M2)3/2]G  (A.  26) 


S4I(xo,  yo),  (x1,  yx)]  -  InKxj  +  V*2  "+fJ>/ 

(xo  +  Vxo  +  y^l  -  G/Vl  +  M2 


(A.  27) 


S5[(xo*  V*  <xl»  yl)]  -  tM^/d  +  M2)372^  - 

(M/(X  +  M2))!  (A.  28) 

where 

G  -  lnK^xJ  +  y2  +  ((xx  +  VjJ/jl  +  M2))/ 

( ^  +  yo  *  <(xo  +  My^/V1  +  M2))!  (A.  29) 

and 

H  -  V*2  +  y2  -  Vx2  +  yj  (A. 30) 

Examining  Figure  1,  it  la  obvious  chat  each  of  the  panel  Integrals 
(2.30)  through  (2.34)  can  be  obtained  by  formulating  the  difference  of 
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tvo  corresponding  .trip  functions.  This  leads  to  the  following 
fundamental  results: 

ri  "  sif(xi*  yi>*  (x3»  y2>]  “  Sl^xl»  *i>. 

<x4»  y2)]  (A.31) 

Tl "  sif(xr  yi>.  <*4.  y2)I  -  s^Cxj,  Yl), 

(*4*  y2)3  (A.  32) 

irhere  1  ranges  from  1  to  5. 
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APPENDIX  B 


SAMPLE  OUTPUT 
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M.PHR  -  5.00 


PRESSURE  RND  SINGULARITY  STRENGTH  DISTRIBUTION 

__ 

_ _ YP  -  «  ISO _ 


LINEAR 


Y/S  -  .038 


x/c 

CRYRE# 

OEL 

CA*» 

cpu 

CPL 

OELCP 

.188 

0.00 

.00215 

.41582 

-.46623 

.36541 

.83164 

.438 

0.00 

.00357 

•193R2 

-.2L074 

.17691 

.38765 

.688 

0.00 

.00R52 

.12030 

-.13151 

’.16910 

.24061 

.938 

0.00 

.00484 

.02656 

-.03436 

.01877 

.05313 

X/C 

cryper 

OEL 

Y/5  -  . 

C&M 

288 

CPU 

CPL 

OELCP 

.188 

0.00 

.00618 

. 38  R96 

-.43403 

.34389 

.77742 

.438 

0.00 

.01164 

.17955 

-.19518 

.16392 

.35910 

.688 

0.00 

.01598 

.11280 

-.12363 

.10197 

.22560 

.938 

0.00 

.01862 

.02556 

-.03343 

.01770 

.05113 

Y7V~= — rmr 


X/C 

CAYRFR 

OEL 

GR" 

CPU 

CPL 

OELCP 

•  1 9  8 

0.00 

“7650  3  T~ 

737  041 

-. 35407“ 

.28675 

.64062 

.438 

0.00 

. 07086 

.14394 

-.15793 

.12994 

.28787 

.658 

6.00 

708090 

.69528 

-.10678 

.08378 

.19057 

.938 

0.00 

.08360 

.01946 

-.02892 

.01000 

.03892 

Y/S'-  778  8' 


X/C 

CRY8ER 

OEL 

CR» 

CPU 

CPL 

OELCP 

.188 

“  5753 

721639 

— 7T3951 

-.38736 

.29166 

.67902 

.438 

0.00 

.17284 

.21860 

-.24551 

.19169 

.43720 

TCTl 

57156 

• 17610 

7F5057- 

-.17154 

7\ 2959 

7361  TT~ 

.938 

0.00 

.18774 

.03844 

-.05524 

.02165 

.07689 

LORO'S  5u Y Yiuy 


YS 

X-CP 

Cho»o 

- Cl- 

C« 

.04 

.28 

2.00 

.48 

.134 

.29 

.28 

2.00 

.45 

.126 

.54 

.27 

2.00 

.37 

.102 

.79 

.34 

2.00 

.45 

.152 
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APPENDIX  C 


PROGRAM  LISTING 


57 


AFIT/GAE/AE/80J-1 


1 

c 

PRC CRAP  NINC<  INPUl,CUTPLT,TAPE5*lNPUT,TAPE6*CnPUTI 

c 

c 

FRCCRAP  NING  CALCULATES  THE  AIRLOADS  Oh  A 

HNC  Of  ARBITRARY  PLANFORP  ANC  CAPBER  IN 

ft 

c  * 
c 

SLeSCMC  FLCW.  A  TRIANGULAR  PANEL  HAVING 

A  BILINEAR  VCRTICITV  DISTRIBUTION  IS  USEO 

c 

c 

TC  PREDICT  THE  PRESSURE  ElELC.  CCPPRESS- 
1PILITV  EFFECTS  ARE  ACCOUNTED  FOR  THRCLGN 

1C 

c 

c 

IN  £  PRANCTL-CLAL'ERT  TRANSFCRPATICN. 

CCPNCN/PLCCKA/XfcC.O  ,V<fcO»2)  ,XCU0I  ,VC<6CI  ,ECtO,2),CV,CX 
CCPPCh/PLCC*E/NS,SSPN,S<15>,NC,CElC),NB,NLI10),B< 1C*2  ),PACH 

CC»NCN/FLCC*C/<<l20,13C)fSG<130>,CBR<fc0),SUM6Cl,ALPHA,NA 
CCPPCN/PlCCKCA  !U5l,ITt;) 

If 

CCPPCN/BLCCNEASPFI5.3) 

REAL  EACH 

c 

"SEAL  IL.n 

2C 

c 

c 

lEClN  INPUT  SECtENCCTCR  GEOpETRIC  OATH 

ALL  CATA  IS  READ  USING  FREE  FCRPAT. 

c 

f E  AC  SPAN  CAT A . 

REACIS,*)  SSPN  ,NS.  (S< I  )  ,I*l,hSl 

IOC 

NRITEIt.lCC) 

FCRPAT<y///2<J>,*SPAN  STATICN?R/26X,*YL0C«.11X,»Y/S*I 

2i 

1C 

“tTTC  I-l.Nl 

NRITEft«110>  I, SII) 

nc 

fc  <FF  mr5>;i3,icy,ri.’n - 

KRIlFCt.UCI  SSPN 

3C 

1?C 

c 

READ  CHORD  DATA* 

RFAT  ( )  NC  ,  It  f  TT,  T* 1,1.  C 1 

NRITEft.1301 

1 3C 

^T**’AT<//2?>,«CFORO  sTa7Icn5*/26X,*VLCC*«11x,*>/C*i 

DC  2C  1*1, NC 

3  * .  ’ 

2C 

INC 

fcWTT E  1 1 , 1  AC  1  TTCTn 

fCRPATI2t>,I?,lCX,Ft.3l 

C 

FlAC  eRfAK  PtlNT~l)ATA. 

RE  AC  IS  ,♦)  NB  ,  IfNLM  ),Bf  1,1)  »B« !  ,2)  ),  I*1»NR  > 

«c 

15C 

NRITF <6,1*01 

FCRNAtf/ /?<!>,«  PREAK  PC  IN  TS* /22  X  ,  *L  E*  ,  ICX  ,  *TEP  ,  10X  ,  *  VLOC  *  I 

3C 

cc  i c  r- itnt 

NRITEIE.IEC)  MI,1),B(I,2),NLII1 

JfcC 

C 

READ  CCNTRCL  PCINT  LCCATICNS  AND  PACH 

4* 

»FAr<S,«l  c>»cp,pXCF 

NRITEU.nO)  PACH,CY,CX 

17C 

“re's  *  A  I  <  /  /  i  7  >  , «  »  A*?H  NUP&Tfi  .♦,Ft.377ITVT»CV-^V,E6.3/31X. 

1«CX  ,Ff .3 ) 

■c 

C 

C 

"  READ  THE  NLPPErTF  LOADiNC  CTSTSTCT? 

FXAPINEC.  THE  ANGLE  CF  ATTACN  ANC  C AP BE R 

c 

c 

- ETrmrcTrcN- fcreach  case  niu  pe  reat 

LATER  IN  TNE  PRCCRAP.  THIS  IS  DCNE  TC 

c 

- FFTSTRVF  S1XR »CE“ SPACED 

RE  AC  <*  ,*  )  NA 

9  W  — 

c 

l>RITEli,«lNA 

CALCULATE  T  ANC  N.  "  IS  THE  NUPBER  OF 

T - CKftNlS'P  f ANTIS  AND  N‘ I  S' THE  Vun'BTT 
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P-NC-  i 
K-KS-1 

*p  -  r  *7 


c 


KPZ*2*KP*P 


CP  SPAM<  I S  E  PANELS 

CALCULATE  PUW6ER  OP  PANELS  HP. _ 

CALCULATE  PUMBEP  CP  UAKNCliK  5  IKGLL AR  IT  IES 


CALL  PESH  CENEPATCR 
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1  SL8CCLTIKE  PE5MKM 


c 

SLFBCUTIfE  TC  CILCULITE  PANEL  CCRfER  PCIK 

t 

Vi  tU  ItKTPCL  PCIKTS. 

CCRfCf/8LCCKA/MfcO«Al  »Y(tC.2)tXC(tO),YC(6C)*E<tO»2l»CV,CX 

"S'  "  "" 

CCrfCf/PLCC«P/fS.S$f>N.Sni>,fC.CUC>,fe,fLU0>,BC10,2>.fACH 

NC«1 

N81-NE-1 

fCl»fC-l 

1C 

CC  <C  V'UKtl 

CC  1C  1-1. fC 

1C 

En,ii«e<K,ii«tin«(B(K.2i-etK,i)) 
flt*2M(Ml  1*  ( B  (K»l, 2  t-e(K«l,  lll-EI  1,11 

IC 1  *  ft  1 K  1 

IC2>M(*«II 

i; 

IC3-!£2^1 

CC  3C  J-1C1.K3 

c?-'  ( it  n-?nm 

C3-SOC11 

CC  Hr-i.sa 

2C 

X  IfC,l)-EU,ll«EU,  2»*«s  CJ  >-03  1/02 

PTTt Ilfrif tnimm  -  .  '  ^ 

xtKc»«i-Eim.n«E»  i*i.2  )♦««*;.  t  >-1:31/02 

YlfC.ll-SUMSSPf 

2! 

Y  ( fT»  2V>$(J«n«SSM 
fC*»C«l 

TT 

1C 

~ rcfrrrfTT - - - 

CCM  I  M  E 

2C 

c 

AC 

“CCMTME 

C JLCULATE  CONTROL  POINTS. 

TC  5C  I  •  1 .  f  F - 

Ycm-n-cv»*ui,ii*cY«vti,2i 

1C 

*c  i»a-CYT*yiT  ,n*  ct*tcy*<xm,a»->it.3»»«ii-cyi* 

1(M!,2>-X(t.m> 

35  ” 

c 

c 

"FrINT  4  T  A  Fl" E- CCT5TS T fMTTE- THE'  PMEL  M3. 

If C  THE  CCRRESFCfCINC  CORKER  AKC  CONTROL 

c 

HINTS. 

►RITE  IE .ICC  1 

ICC 

"TCfff  Mt//y/tX.*PAfEL*»T>.*Xl*.Bli.«X2*fflX,*X3*,?y;*XA*»"^y7BVl*T 
10X  t«>2«.9X,*XC«.8X,PYC*1 

rCTt_T^Tf>P - 

minu.icil  l.»Utn.mi.2l,)M1.31.«U.A»«TM,H  .YU.2I.XCU  >. 


- jYcrn - 

101  t* ,12, f (4X.F7.2) » 

r* - sr  rrwmnr* - 
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1 


1C 


1? 


c. 

c 

c 


c. 

c 

c 


1C 


SLBfiCt  TIKE  IM  (>liX2i>}iliMl«f2) 

_ U.I5_SL'BRCITIKE_  EVALUATES  TFE_PAKEL  Hllfc. 

cfals.  It  als'o  flkcticks'as  theexecltive 

_ CCMBQL  RCLTIKE  FOB  SueRCUTIKE  STRIP. _ 

COPPCK /PLCCKCA It  15  I »I TIE l 

CCPtCK/PLCC*E/SPF(5<3> _ _ _ 

REAL  II, IT 

_ CALCULATE  S  TRIP  FtlKC  T  IONS  F  OR  THE  PAKEL 

LEACIKC  EOCE  KITH  CORKER  POIKTS  m»Yll 
AKC  1X3, Y2). 

I-l 

CALL  STPIPQ1,Y1,X3,T?,1) _ _ _ 

CALCLLATE  STRIP  FUKCTICNS  FOR  THE  PAKEL 

_ MU  OIACCKAL  KITH  CCBKEB  PCIKTS  <X1,Y1I 

AKC04.Y2I. 

J-2 _ 

CALL  STRIP01,Y1,XA,Y2,I» 

CALCULATE  STRIP  FUKCTIOKS  FOR  THE  PAKEL 
T*  A  IL IKC  ECCE  WITH  CCRKER  'POINTS  IX2,Y1I 
AFC  IX*,Y2I. 

Tv5 - 

CALL  STRIPf>2,Yl,*4,Y2,tl 

- pyatia  tn  AKfinKT  exralst - 

CC  1C  1-1,5 

t  cnnim  TTir^miTn - 

nm-spF(i,2i-sPF(i*3> 

CCN'tTmJe 

SETLSf 


2C 


C 

c 

c 
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1C 


rr 


2C 


Tf 


3C 


*W 


j 

i 


AC 


; 


i 

i 


SLP'CLTINf  STMP<XO.YC,>ltYl,I» 

C _ SLPRCUTIKE_TC_jCALCUl»TE  STRIP  FUPCTICN5. _ 

C  TPERE  ARE'  15  STRIP  FLKCTfCK  CALCLLATIOTS 

C  PER  CU»CRIL*TERAL  PAPEL. 

CCPPCP/BLtCPE/SFF<5,3»  ' 

RE  * l  PC 

T  isnTPff  fl-pCTTcp  ecu* t Tops  for  an~Tk- 

C  PIMTE  SLOPE  IF  Pl-XC. 

tFI*P3f»I-XC».LT.ie-ei  CC  TC  ICC 
PC-  Ol-YCI/tXl-XCT 
BC-YC-XC4PC 

CC»1-SC«T(>C«*24Y0**2» 

OLP2-SC*r«>I*«i«YlP»2l 

CUP3-SCRTM«PC«*2> 

Cu  p  * • t  CL  P i • ( x  c «  p  c  * y c  Yit up 4 1 7TftP7TT x TV pcTv i > / c  L  p3*T 

CLP4-/ESICLP4  I 

ctvv-  h  c  c  <  c  l  pvt 

CLPi-iPStCUPlOCI 

cl  p*v*nrcmrm 

CLP<-*PSCCLP2«>1» 

Cl, e-iicciciPtl  ■ 

SPFlI^n  — PC«CLPA/0LP3 

S>R  t<Tf»*YC«CLP'5-Yi«CLPe~PO»TcUPT^CUKI)/DL'PT*T2-RT‘tTPA/OL,P3*P3 
5PF  13  .II-»C*«<*ICLP?-Cu»n/CLP3»*2-PC*P0»Cl;H<,/CU13P«3 
SP'F'f  a,  1  i- CUTT/ C  t  p  3  ♦  t  L  p  e  -  C  Lr"  5 

SPM! ,11— eC«PC**2*nLP*/CLP3«*3«PC*IDLPI-DUP2>/DUP3«*2 

- RlTul^ - 

ICC  CfMIPLE 

- Clpi-sc cTTvr*  r?v>c*T?i - 

CLP2-SCRTOC*«2»YX*«2» 

- clt^*  irsrrrrPT*  yttt7  Trcrirrm - 

SPFtl,l»  — UCCICUP3I 

CvRirrFSut^ioc)  “ -  ~ — 

OLP3-»ltC«CLP3I 
OlT*.  *1  ST  CTT»  2 « rn 
CLP4-IICCICLP*  I 

Trt  li,  II-v£«ClP3-Yl+CUPA  "  ■  ' 

SPFI3.H-CL'P2-tOPI 

- sPFivm-tov-ttRj - 

spft«.ii— *c*SRFu.n 

irrmfR  ■ — ~ — ~~ 

epc 
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w 


WWW 
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1 

C 

SIBSCITINE  1MFTIPP2J 

SUBROUTINE  TO  INVERT  THE  AERCOVNAPIC  IP- 

c 

c 

flUEPCE  CCEFFICIEKT  ARRAY  USING  GAUSSIAN 
EUP1NATICN. 

• 

CCHPCMPtCCRCAA (130*1301 «SC( 130 ) ,CBR(6C 1 «SUP( (Cl  *  ALPHA, NA 

CC  2C  I-1.KF2 

pivcwu.n 

Ml.tlM. 

1C 

1C 

CC  IC  L-l. NB2 

AH  ,H-A  (  1  .LI/FIVCT 

rc  2 t  P-1.PP2 

TFI'.EC.I)  CC  1C  20 

ii-nr,n 

a«p,ti«c.c 

is 

15 

CC  15  1-1.PP2 

AIP,t»-AIP,t»-AU.t)ATT 

2C 

CCN11M.E 

RETURN 

ENC 
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lEfCUTlKE  4EFt(KP«f»K,KP2> 


IMS  SUBROUTINE  CUCUL4TES  THE  4EBG0YNI- 
P I ?  INFLUENCE  COEFFICIENTS  ASSUPflNC  * 

\C  RTIC  ITY  CISTRieuTICK  Of  <-eX*ET«F)l« 
4«BX»~CYt  J.  THERF~ *RE  PINT!  UkKNCRHS 
If  TER  THE  FCUNCIRY  CCNOITICNS  H* YE  BEEF 

fortr: 


CCfftf/PLCCKI/Y  160,4)  ,Y<6C,2),XC(60) *YCC6C) ,Elt0.2) ,CY,CX 

ccpp CN/ertcxe/fs  .sspn.susi'.nc.cuo.nb.nluo)  ,Mio  .2)  ,p*ch~ 

CCPFCF/PLCCKC/I(130,130),SC(130)*CBR<6C),SUP(6C>*4LPH4*N4 

TC  7>  'C  ) ')  ETT'f  x  C7 1  LI  STVl  TT51 - 

REU  RICH, II, IT, KL, XT 

"TKTFin! — CTTCTYD  370  ^7017077070  4 - 

PI-3.1415S2 


rc-c 


2ER0  OUT  THE  4ERCCYN4PIC  1 

NFLUEKCE 

ITTon 

CC  IC  I-1.NP2 


fCINT  J  CUE  TO  P4NEL  I. 


DC  37C  I-l.RP 


IPRLY  *  IINE4R  TR4NSF0RK4TI0K  TO  THE  CC- 

Tnrrrrn  rx  mximvt  sr  tf  r~vrr  or  i  m - 

IS  IXCJ.YCJ). 


AFIT/GAE/AE/ 80 J-l 


C 

CELT*? 

C2-! 

cc 

c 

CELTA3 

C3-I-HP 

c 

CE11I4 

CA-t«» 

t« 

c 

«rp»l 

C1-NP«P«I 

c 

(«PP*2 

C2-NF«P«I.l 

c 

(IPP13 

C3-KP«2*P.l 

If 

c 

OPPI4 

C*»NF«2*p«  I«1 

c 

c 

If  ELY  THE  PRANCTL-CLAUERT  TR ANSFCRPATICN 

7! 

c 

XI  ■ X 1 /If CH 

X2-X2/AK CH 

xs-xi/apch 

c 

XA-M/APCH 

{ ILCULITE  LEADING  EDCE  flKEL  CHECK  PARAP. 

c 

ACTICE  THE  USE  CF  INTEGER  ARITHPETIC. 
Nt«I-II/P>AP-I 

c 

CILCULf TT- 3b* IL1KG- iDCE  PANEL  CHECK  PABIP. 

NS-1-II/P1*P 

*« 

t 

c 

A 

THIS  LCCP  CCPPtTES  THE  INDUCED  VELOCITY 

c 

c 

AT  CONTBCL  PCINT  J  OLE  TC  PANEL  I  IK-11 

ANT  THE  IPACE  PANEL  CF  I  |K-2>. 

c 

CC  3«C  K-1.2 

ttP*«-l.»«*» 

Y1-YII,11«YCJ«CVP 

c 

Y2-Yll,2»4YCj«CtP  ’ 

f! 

c 

c 

EVALUATE  ThTTEN  PANEL  INTEGRALS. 

c 

CALL  IHT(Xl,12»X3.X4.Vl.V2> 

Kl-«?MYl-Y2»«ILm-IXl-X3MILHI-IVl*X3-Xl*Y2»*IL«OJ/IYl-Y2I  . 

«t-i?*iyi-y«  »«nm-«>2-*o*nm-m*x4-x?*v2»*iTnn/m-Y2) 

iff 

c 

CHECK  IF  PANEL  I  IS  A  KINC  TIP  PANEL. 

IF  1 I.CT.f XP-P1 1  CC  TC  2C0 

c 

c 

MNEL  ECIATICNS  IN  THIS  SECTION  CF  THE 

ic« 

c 

c 

ftCiTlE  FCR  BCCT-ThCBC  ANC  INTER 10R  PAXElT: 

c 

IFCKJ.KE.PCI  CC  TC  2<J 

AFflY  FLOP  TANCENCY  CONDITION  AT  THE  L.E. 

Y3-VU«P»2l«YCJ*CtP 

XS-IMHP.3IOC.JI/APCH 

c 


X  ‘  ■' 

_ aij,C3i-ai  j.c3t«HLisi-Ti«iLMn»ix«..A3>/m-T?i/iY3-Y?>/«m 


lie 


AFIT/GAE/AZ/80J-1 


1 i; 


12C 


12! 


11C 


!«C 


TXT 


1!C 


TTT 


IfC 


ITT 


2£  CCM!»U 

C _ _  _  _ 

t  <  j ,?  1 • •  *  i  j  .1 1  i « « »/•  ii  <  *  T-  It  i*T»/  itj -i?  i  /*/n 

c _  _ _ 

c 

Ic  mffVCF 

c _ _  _ _ _ _ _ 

•li.mxu.d  i«  tt  *r*>2  tVnii  i-i»2-mi*v2*iL  <i>-m-Y2i 

_ I»HI2MU->|  l_/«»2-»il/iY|-Y2I/42M _ 

C 

_ >Uit2i»n  j  .c;  hhmiiio-  t  tun/m-tn/wn _ 


c 


c 


c 


T 


c 


2CC 

x — 

C 

T - 

C 


c 

T 

c 

T 


22C 
T - 

^23C 

T - 


e 


Min«».>|»»)«  lHtt«tY2-YIf  lH2»«Qil-*4f  «H 

T777 <  -  y 1  i  / « *  ?-*« » / »  /  rf~ 

<  uTT '  >•  ti  j  ,C*Th  lniMffTf  i«  »f/l>l-V2774/>  I 
"TTTRm .>T TT T  Tf  Jit - 


iTTTTTTol  j.et  i«  H  v ? •  > l- 1 l  *11  *  ii«4  V i-«ri » *  Hi 2  TOT 

l/IYl-Ytt/t»<- 


>H/4/H 


J« j,CM*<n0«v;-Yl*>]>*lLll>-IX3->O*Yl*M U»-tY2-VlJ 

T*  IL  1 2  I « I *1- M  M  Kl» I *1-X  4  l*K  TTm  Y2” YllTl  X1-X4  I  / 4/1*  I 

-rrrnn: - 

CCkllFUE 


MMt  fCtillCM  tv  TvIS  SECT IC*  Cf  TVF 
- rrcFTsr f r *“TTt^cHcrD  mkfls - 

trnr:LETn>  »  Tt  ?7c - 

/FPLY  FLOW  T4NCEVCY  COMHTICfc  t T  1HE  L.E. 

- M  TF  -  C*fX*3~C.“TmTT$~tlLTiT~^_0^ - 

MICE  L.E. 


*  t  J  <C1  *•*»  J.Cll«<Y2*Kt«»-TH5n*l>1-»n/IYI-Y2»/tY2-Yn/V/fI 

Cc  It  tie 

CCMIMiE 


j.en*(Y2«iu«>-tLt;>i/iYi-Y2t/4/Pi 


•<  j,cn*tiLm-Yi«iit4t  i/«n-Y2t  24/pi 


CCNYITVF 


J.Ct)«lY2*IT(«l-n<H  >/<Yl-Y2t/4/Pl 


M  JtC4».««  J,CO*tITt5)-Yt«n<41  l/«Yl-Y2>/4/f  t 


<f  J.Cn*ft»44Yl-Y24>2l4ITm-tX2-*n4Y2*ILtn-tYI-Y2» 

1«IT  f  2  l«t>7-VTHKL«»X2-X*  |4KT » / f> 2 -XTT/fYT^ T2T7 *7 P  l - 
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c 

WIY1-Y2I/02-P1I/4/PI 

175 

sec 

37C 

CCM1MJE 

fCMIUf 

•  4CC 

c 

CONTINUE 

c 

c 

IMS  SECTION  OF  SUERCUTINE  AERC  FORMULATES 
TPE  NP  CELTA  SPECIAL  CONTINUITY  CCNDITICNS 

1<C 

t 

c 

AT  TNE  LCNER  RICM  HAND  PANEL  CORNER. 

CC  tCC  I-1«NP 

INP-HNP 

ies 

xi*>ii»ii#apcn 

X2-MI.2I/APCP 

>3011  *3  1 AAPCN 

X40IT,4I/APCF 

V  10(1.11 

Y20II.21 

ISC 

c 

C3-1-1 

C4-1 

Cl"NP«l 

155 

c 

niCuLAY?  It  ApC  TTTnEcf  parameters. 

NE-I-II/MlAP-l 

c 

N5.|-;1/M)«p 

CNECX  IF  PANEL  I  IS  A  RING  TIP  PANEL. 

2CC 

c 

Tru.n.  INP-Mi»  CTTTO  5CC 

c 

c 

Panel  CONTlNuITf"EouA"TropS_Tw  This  stCTTBN- 
CF  THE  CODE  ARE  FOR  INTERIOR  AND  ROOT 

205 

c 

c 

'  CRHUfFINetT: 

TIT 


410 


420 


215 


AIINP.C3I  —  IX5-X3IMY3-Y2I 
CC  TC  420 

CONTINUE _ 

■7TYNP7153I.-1. 

CONTINUE 

TTTkp. 041-1. 

aiinp,c3i«aiinp,g3)-ixi-X4>/iy2-yii 

TrjTTTfrrn  >  so*  ft  ivz-yi  i 

If  IPS. EO. PCI  CO  TO  590 _ 

*n4P.c4»— Txa^Fimr^fn 
CO  TO  590 


— Tf(i - 

C 

C 

PANFL 

CONTINUITY 

EOUATIONS  IN  THIS  SECTION 

C 

c 

CF 

THE 

CODE  ARC 

Po’r  TIP  c'hcTro  panels. 

225 

5w> 

“EBnTTnuE 

IF  f  N8.E0. PC  1 

o 

TO 

510 

510 

AUPP.C3I  — 1 
CONTINUE 

• 

- AC INP .04 l • 1 . 
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A<INP(C1I-<X3-X4)/IY2-V1> 

3<?0  CONTIMIt _ 

600  CONTINUE 


THIS  SECTION  OF  SUBROUTINE  AERO  FORMULATES 
P  ADDITIONAL  DELTA  CONTINUITY  CONDITIONS 


CN  THE  ROOT  CHORD  PANELS. 


DO  TCC  1-1. P 
INP«2*NP« I 


J-I 

Xl-XU.il/APCH 


X2-XU.2I/ABCH 

X1-XU.3I/APCH 


X4-XI J.AI/APCH 
Yl-YU.ll 


C2-NP-PO-1 


C*-NP*2*H-»J-1 

CALCULATE  TE  ANO  LE  CHECK  PARAMETERS. 


Nfl-J-I 

N9- J-I J/-IPP 

-jrfKe  .“nE;>'CI  "CC'TO  710 - 

AIINP.C1I— IX3-X1J/IY2-YII 

“Cff~TD~72I5 - 

CONTINUE 

-rrnsPTo  i»--i. - - 

CONTINUE 

TT I N  PT02T=TT - 

AUN'P,ClI-AUNP,Cll*IX2-XAI/CYl-Y2l 
I F  t  N  ^  •  E  0  .MCI  C  tr  XWYOO 
At  IMP.  CAT—  IX2-XH/IY1-T2I 

~nTN>Tc2r-(kA-ia»/<Vl-Y2> - 

CONTINUE 


AFIT/GAE/AE/80J-1 


SUBROUTINE  LOACS<NP,M,n,NP2I 

_ SUBROUTINE  TO  CALCULATEJPRESSURE  D1STRI- _ 

PUT  IONS  ANO'aERCDYNAmIC  COEFFICIENTS. 

_CCNM0N /BLOCK  A/XCfcO,jO jYi60.2».XCI60I.YCC60I.E  f  60, 2 1_.CY  .CX _ 

CONMON/RLPCKB/NS.SSPN.Sf  15 1 ,C (10) ,NS ,NL f 101 ,8< 10 ,2 ) .MACH 
_COM  «ON/RLO.CJ<C  /  A 1 130 .130)  .SCI  1 30 >  •  CBR I  b_0 >  •  SUM  160 J  .AL PHA.NA 
DIMENSION  0EL(60).GAM(60).CPU(60) .CPL I 60I.CLL (60 I ,CMf 60 I 

INTEGER  Dl.D2.03.0AtGl.C2.G3.CA _ 

REAL  «ACH 

C- 1. A _ 

A«CH»S0RTIl-MACH»*2 I 

PC-0 _ 

SET  CAMBER  SLOPES  EQUAL  TO  7ER0. 

on  10  l-i. np _ 

CRR  f  1 1"“0  . 

CONTINUE _ _ _ 

READ  NA  SETS  OF  ANGLE  OF  ATTACK  AN0 

_ CA»BER  SLOPE  PI STRIBUT10N  DATA. _ 

00  800  L* 1 »  N A 

READ(S.P)  alpha, nccg.npf.cy.cx _ 

' _ NCCC_IS  THE  CAMBER.  CHANCE  PARAMETER^ _ 

ENTER  0  TO  REAO'a  NEM  CAMBER  SLOPE  DISTRI¬ 
BUTION  0«  ENTER  l  TO  RETAIN  THE  PREVIOUS 
DISTRIBUTION. * 

npf  is“t_ bTT ressure  >*tion~ par  a  me  TerT 

_ ENTER  0  TO  USE  THt  EXACT  ISENTROPJ[C _ 

EXPRESSION  OR"  ENTER  l  TO  USE  THE* LINEAR- 
_ 1 2CD  FORM. _ 

COMPUTE  PRESSURES  AT  PANEL  LOCATION  CX 
tNftY.  NOTE"-  PRESSURES  MAY  bETOmpuYFO 
_ AT  POINTS  OTHER  THAN  THE  CONTROL  POINTS. _ 

CALCULATE  PRESS  IRE  EVALLATICN  PCIPTS. 

"EC  If  r-  I  .  PI 

3cni*<i-c>t<Y<i.n«CY«Mi.;i 

>(« )1«C3«> 1 1,31«  Il-CYIO  «|,1  l«7>T(0«  |>  1 1 1 A  i  -  TTTTT7T«TT*7y7« 

i  t)iitii->iii  ini _ 

I  (I  M  (  (  .  ft  ,  I  1  ((  ic"7c 

RfIMf.O  KMHI.fl.Ml 
Cf»n»ll  ' 

PtFHtMf  lIH/Mlft  Ft  F  F  CF  TH  FlCV 

rotifO  rrociFi  conno.  ftvni 

•F  I  F  FI  SI  MS  IF  E  FlFPUIlEt  7  CCPFIFEI  CF 

ufinni  rmifii  r. 

C(  :(  I  •  1  ,M _ 

St  M  1 1  •  1  CTTTTToTffTTTS 

cmiui _ _ 

uminin  fiF  TTTn  ey  /  H  ce_imn  TFI 
tFFFEFF  SlFCtUFIIY  SIFEFCIFS. 

It  f  t  J •  1  . M  i 


tt  •( 

hf»«, 
t  '(  ‘  c 
root 

Tnrnr 


!• 1 »PF 

MMJ,n<SLMII 

n - 
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«( 

cmmt 

<( 

f 

IH1M. 

C  - 
( 

«wm  itit  fcf  ccpfimc  m  simiAFin 
nmm  n  im  tmm  mu. 

e* 

C 

CC  <CC  I-l.AP 

c 

c 

FETPIEYE  CCMRCL  PC  I  Ml  AKC  PAKEl  CORNER 

EC  IMS •  APPLY  FRAKCTl-CLAUERT  TRANSFCRP. 

xi*xct  i  i/amch 
ri-vcm 

70 

xi-xiuiwapch 

X2-XI! .21/APCH 

X3*XtI.3)/AKCH 

X4-XI  I *4) /APCH 

75 

Yl-Ytl.l) 

V2-YII.2) 

c 

'  ‘  CALCULATE  THE  8  SINGULARITY  NUNBEfcS. 

01-1-1 

02-1 

03* I-1*H 

BO 

0 A* I  ♦  H 

Cl«NP»tt«I 

C2*NP ♦  •*♦!♦! 

85 

c 

G4*  NP*2  I  ♦! 

CALCULATE  TE  AND  LE  CHECK  PARAMETERS. 

NS*  t-l  !7*n •  «ST 
N9*t-II/N»*N 
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GA-4-0. 

DEL  1-0* 

95 

DEX2-0. 
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c 
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c 
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110  CONTINUE 
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OEL3-1GID3I 
tto — cwttnof — 
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CA93.SG<G3l 

IFIN9.F0.MCI  GO  TO  900 

GAN?-SG(G2l 

m _ 

GAN9-SG(C9t 
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c 

CONTINUE 

c 

c 

SINGULARITIES  DETERMINED  IN  THIS  SECTION 

ARE  FOR  TIP  CHORD  PANELS. 

125 

c 
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DFLI-IX1~X1I/IV2-Y1 )*SC(G1> 
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310 

CONTINUE 

DELl-SCfDl) 

320 

DEL3-SGI03I 

CONTINUE 

0EL2*$C(0?1 

0EL‘-SGf09l 
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GAMl-SGICll 

IFIN9.E0.NCI  CO  TO  900 
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CONTINUE 
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c 
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c 
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c 

BL-IGAN3-CA«9»/(X3-X9» 
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1  /  I  X  2— XI 1 / 1 Y 1— V2 1 
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c 
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C  START  OUTPUT  SEQUENCE 

230 _ C _ _ 

W«TTE(6,610>  ALPHA 

_ 610 _ F OR  MA TJ^ / 1 tJ2* X . PALPHA  -♦,F6»2> _ _ 

VRITEI6,620) 

_ •  620  FOR  M  ATI  /  / 1  2  X_,  ♦  PRESSURE  AND  SINGULARITY  STRENGTH  DISTRIBUTION*) 

235  WR  I TE 1 6*625 1  CX.CY 

_ 625  F0RMATI/30X.PXP  ■ ♦ , F 6 . 3 / 30X , *YP  -P.F6.31 _ 

lFINPF.EO.il  GO  TO  627 

_ WRITEI6.626J _ 

626  F0RMATI/33X.PEXACT*! 

2M0 _ _ CO  _T0_629 _ ; _ 

627  CONTINUE 

_ WR_I  TF  1 6. 62  8J _ _ _ 

628  FORMAT I/32X»*LINEARP) 

_ 629  CONTINUE _ _ 

2M5  Nl-1 

_ 00  760  J-1,N _ _ 

6u“»su»«cv*csi  n-si  jn 

WR  T  TE 16 1 630 1  OUM 

630  F0R»AT(//30X,'»Y/S  »*.F6.3/  6X ,* X/C*. A X. *C ANBE R*, 5 X, *DEL *«6X ,• GA «♦ , 

16X,*CPU«,6X,*C?L*.6X,*DELCP*I 
00  750  I-l.M 

OU*l-Ctn«CX*(C(I*ll-CIMl  _  _  •_ 

~t)UM~2»CPL  I NII-CPUI Nl I 
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750  CONTINUE 

_ 760  CONTINUE _ 
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c  _  _ 
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00_  790_x -lfNBl _ 

Nl-NLIKI 
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oum-si jTTCypisi j«ii-si j >T 

0UM).|0U«-S(N1I)/IS(n3I-SIN11I 

Bum  2-  *P*K,2  >-BIK,ll  IPIl-DUMl  I  ♦Bum*  fB  ( K»l,2I-6(K«l«  II I 

DU  M3  -  0  . _ _ _ _ _ 

DUM5-0. 

_ 00  770  I-l.M _ _ _ 

NA-NsTi 

280  0UM3-0UM3*CLLIN5)  _ _ _ 
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770  CONTINUE  _  _ 

0UM5-0UM570UM3/DUM2 

OUM3-nUM3/OUM2  _  _  _ 

"285  '  OUMA- OUMA/1 DU M2**2) 
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250 
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WRITE (6.7751  DU*' , DU 65  . 0U*2 ,  OU M3  ,0U*6 

_ 775 _ fO  R  *A  T(l<,X.f5.2.3X<F5.2.2X.F 7.  2 .2 X  .F 7. 2 , 2X.F 6 „ 3> 

780  CONTINUE 

_ J790  CONTINUE _ 

290  800  CONTINUE 

_ STOP _ 

END 
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